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Abstract 

We investigate theoretical interpretations of the 1.4 GeV J PC exotic resonance 
reported by the E852 collaboration. It is argued that interpretation in terms of a 
hybrid meson is untenable. A K-matrix analysis shows that the 1.4 GeV enhancement 
in the E852 rjn data can be understood as an interference of a non-resonant Deck- 
type background and a resonance at 1.6 GeV. A final state rescattering calculation 
shows that the 1.6 GeV hybrid has a r/ir width which is bounded above by 57 ± 14 
MeV. 
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1 Introduction 


Evidence for a J FC = 1 h isovector resonance p(1405) at 1.4 GeV in the reaction 7T~p —> 
pir~p has been published recently by the E852 collaboration at BNL |]J. The mass and 
width quoted are 1370 ± 161™ MeV and 385 ± 40 ^io 5 respectively. These conclusions are 
strengthened by the claim of the Crystal Barrel collaboration that there is evidence for 
the same resonance in pp annihilation with a mass of 1400 ± 20 ± 20 MeV and a width of 
310 ± 50±|8 MeV ||, consistent with E852. However, the Crystal Barrel state is not seen 
as a peak in the pir mass distribution, but is deduced from interference in the Dalitz plot. 
Since the J PC of this state is “exotic”, i.e. it implies that it is not a conventional meson, 
considerable excitement has been generated, particularly because the properties of the state 
appear to be in conflict with theoretical expectations. The resonance is reported in natural 
parity exchange in the E852 experiment, and no statement can currently be made about 
its production in unnatural parity exchange. 

In addition there are two independent indications of a more massive isovector J PC = 
1 _+ exotic resonance p(1600) in ir~N —> ir + n~ir~N. The E852 collaboration recently 
reported evidence for a resonance at 1593 ± 8^7 MeV with a width of 168±20t}2° MeV [[|. 
These parameters are consistent with the preliminary claim by the VES collaboration of 
a resonance at 1.62 ± 0.02 GeV with a width of 0.24 ± 0.05 GeV (4j. In both cases a 
partial wave analysis was performed, and the decay mode p°ii~ was observed. There is also 
evidence for p(1600) in r/V peaking at 1.6 GeV 0. It has been argued that the pn , p n and 
pir couplings of this state qualitatively support the hypothesis that it is a hybrid meson, 
although other interpretations cannot be entirely eliminated 0 . 

Recent flux-tube and other model estimates 0 and lattice gauge theory calculations [§] for 
the lightest 1 h hybrid support a mass substantially higher than 1.4 GeV and often above 
1.6 GeV ||. Further, on quite general grounds, it can be shown that an pit decay of 1 -+ 
hybrids is unlikely [0 . There is thus an apparent conflict between experimental observation 
and theoretical expectation as far as the 1.4 GeV peak is concerned. 

The purpose of the present paper is to propose a resolution of this apparent conflict. Two 
possible hypotheses are considered. 

1. The two states are indeed separate resonances and are hybrid mesons: the lower one 
the ground state and the upper one an excited state. We perform calculations in the 
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flux-tube model of Isgur and Paton |1(| to demonstrate that both on mass and decay 
grounds, this hypothesis is implausible. 


2. We suggest a mechanism whereby an appropriate t/tt decay of a hybrid meson can be 
generated and argue that there is only one J PC = 1 _+ isovector exotic, the lower- 
mass signal in the E852 experiment being an artefact of the production dynamics. We 
demonstrate explicitly that is possible to understand the 1.4 GeV peak observed in r/Ti 
as a consequence of a 1.6 GeV resonance interfering with a non-resonant Deck-type 
background with an appropriate relative phase. We do not propose that there should 
necessarily be a peak at 1.4 GeV; but that if experiment unambiguously confirms a 
peak at 1.4 GeV, it can be understood as a 1.6 GeV resonance interfering with a 
non-resonant background. 


2 Hypothesis I: Two hybrid mesons close in mass 


The simplest explanation for the experimental report of two peaks at two different masses, 
is that they are indeed separate Breit-Wigner resonances. 


The most conservative assumption is that these are then both hybrid mesons. Other less 
likely hypotheses, such as glueball, four-quark and molecular interpretations, are discussed 
in ref. 


In the hybrid scenario, the 1.4 GeV resonance would naturally be assumed to be the ground 
state hybrid and the 1.6 GeV resonance an excited hybrid. A numerical calculation in the 
flux-tube model indicates that the orbitally excited D-wave hybrid is the lowest excitation 
above the P-wave ground state, with a mass difference of 400 MeV for light quarks [|7j. 
The same model predicts cc D-wave hybrids to be 270 MeV heavier than the ground state 
hybrid 0, in good agreement with the result of 230 MeV found in adiabatic-limit lattice- 
gauge theory simulations [§ and similar results in NRQCD lattice simulations [|lTJ . Also, 
the lattice-gauge calculations find that the next highest levels in the cc sector are the 
radially excited P-wave hybrids [|TT] , which are 400 MeV heavier than the ground state [Q. 
A mass difference of 400 MeV for the light-quark hybrids is clearly inconsistent with the 
experimental claim of resonances at 1.4 and 1.6 GeV. The absolute mass scale predicted by 


3 











Tabic 1: Decay widths of a ground state at 1.4 GeV in the flux-tube model. The 
/ 2 7 t, pui, K*K, pn and pn modes are all substantially below 1 MeV. The conventions 
and parameters are those of ref. |H|, except for the following changes for the bin and fin 
modes. Here we use a radial dependence of the hybrid ~ r, which produces widths ~ 5% 


different from ref. |L3|. More importantly, we take into account the fact that the b\ and 
f\ have finite widths, and we assume that they decay predominantly to un and a 0 (980)7r 
respectively. 


Decay Mode 

Partial Wave 

Width (MeV) 

bin 

S,D 

96 

/w 

S,D 

13 

pn 

P 

4 

77(1295)tt 

P 

< 1 


theory does notj] support a ground state hybrid at 1.4 GeV, as discussed in ref. Thus 
there are two arguments on mass grounds for discarding this hypothesis. 

Further, from the viewpoint of decays, it is qualitatively hard to explain why the lower- 
mass 1 _+ hybrids should be seen only in r\n. This is because relativistic symmetrization 
selection rules suppress the pn decay of any 1~ + hybrid in QCD in the absense of final 
state interactions ||. Within the flux-tube model and constituent-gluon models there is a 
selection rule which suppresses decays of ground state hybrids to two S-wave mesons Jl3]. [14 


This selection rule requires only the standard assumptions of non-relativistically moving 
quarks and spin 1 pair creation in a connected decay topology ]TB[. In addition, for 1 h 
hybrids the selection rule is only operative when the non-relativistic spin of the QQ is 1. 
The lowest orbitally excited hybrid in the flux-tube model has QQ in spin 1, and hence 
obeys the selection rule. The ground and lowest excited hybrids have hence got similar 
overall decay structure. 

Flux-tube model predictions for the decay of a 1.4 GeV hybrid are given in Table 1. We 
note that the total predicted width of ~ 110 MeV is much smaller than the observed value. 
The calculations show that we expect an appreciably larger pn width than pn width for 


the ground state hybrid. This is confirmed by QCD sum rule calculations |16|. It then 


1 Except for a QCD sum rule prediction of ~ 1.5 GeV 
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Tabic 2: Decay widths of an orbitally excited hybrid at 1.6 GeV to P + S'-wave states 
in the flux-tube model in MeV. The conventions and parameters are those of ref. |13]. 
The derivation of the widths is discussed in Appendix [A|. The inverse radius of the hybrid 


/ 3p = 0.27 GeV is taken to be the same as that of the ground state hybrid 


We also 


quote an error based on taking /3p = 0.23 GeV. Widths of a 2 GeV orbitally excited hybrid 
is also indicated for comparison. 


Decay 

Partial 

Excited Hybrid Mass 

Mode 

Wave 

1.6 GeV 

2.0 GeV 

biTT 

S 

118 (-22) 

10 (-5) 


D 

.1 

.8 

/w 

S 

30 (-5) 

4 (-1) 


D 

.05 

12 (-1) 

fl* 

D 

.08 

2 

a\T) 

S 

- 

11 (-3) 


D 

- 

2 

a 2 r] 

S 

- 

.07 

Ah(1270) A' 

S 

- 

105 (-27) 


D 

- 

.5 

Ah (1400) A' 

S 

- 

.05 


D 

- 

.3 

AT* (1430) AT 

D 

- 

.03 


becomes difficult to understand how there can be almost no presence of 1 _+ wave in the 
pTT experimental data at 1.5 GeV ||, where there should be significant presence due to the 
~ 400 MeV width of the E852 1.4 GeV state. This calls into question the interpretation of 
the 1.4 GeV state as a ground state hybrid. When final state interactions are taken into 
account (a point on which we elaborate below), we expect a larger ijtt width, which may 
invalidate the preceding arguments. We shall hence proceed with the hypothesis that the 
1.4 GeV state is the ground state hybrid and the 1.6 GeV state the orbitally excited hybrid. 

According to the flux-tube model calculations in Table 2, the orbitally excited hybrid at 
1.6 GeV has a somewhat larger total width than the the ground state hybrid at 1.4 GeV. 
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This is a strong theoretical statement as generally nodes in orbital wave functions tend 
to suppress specific partial widths relative to the ground state. Note that P-wave modes 
like 77(1295)7r, K*K, pu should all be strongeif] for a 1.6 GeV state than a 1.4 GeV state, 
simply due to phase space. Thus there is a further problem in understanding why the 1.4 
GeV state should have a larger experimental width than the 1.6 GeV state. 

So on a multiplicity of grounds we are forced to conclude that the hypothesis that the 1.4 
and 1.6 GeV states are both hybrid mesons is theoretically untenable. 


3 Hypothesis II: A single hybrid meson at 1.6 GeV 


The current experimental data on the 1.6 GeV state is consistent with mass predictions 
and decay calculations for a hybrid meson 0. O. This then leaves open the interpretation 


of the structure at 1.4 GeV. 


There are two basic problems to be solved. Firstly it is necessary to find a mechanism 
which can generate a suitable ijir width for the hybrid. Then having established that, it is 
necessary to provide a mechanism to produce a peak in the cross section which is some way 
below the real resonance position. 


We first show that a sizable 777 r width for a hybrid resonance can be generated by final- 
state interactions. For this we use a doorway calculation, the procedures for which are well 
established ||T 8 |]. We use the simplest approach to provide an upper limit. 

The 777 r peak in the E852 data spans the /wr and 617 r thresholds, so we propose a Deck-type 
model [0 as a source of a non-resonant 7777 background. We then show that, within the 
K-matrix formalism, interference between this background and a resonance at 1.6 GeV can 
account for the E852 r/n data. The width used for the decay of the 1.6 GeV hybrid to 7777 
is comfortably below the upper limit established in the doorway calculation. 


2 In constituent gluon models, the lowest-lying excited hybrid is expected to have QQ spin 0, so that 
decays to S-wave mesons are not suppressed Bill- The lowest-lying excited hybrid would then be very 
wide indeed. 
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Figure 1: Decay of p to pn via final state interactions. 
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3.1 t]tt doorway width of a 1.6 GeV state 


Although the pn width of a hybrid is suppressed by symmetrization selection rules 0 which 
operate on the quark level and have been estimated in QCD sum rules to be tiny (~ 0.3 
MeV) jn|, long distance contributions to this width are possible. We shall show that these 
can be very much larger than the widths obtained without final state interactions. 


The procedure we adopt is that of a doorway calculation with on-shell mesons [18|], which 
provides an upper limit. An essential ingredient is the presence of an allowed dominant 
decay which can couple strongly to the channel of interest. In the flux-tube model b\n is 
such a dominant decay, and it is strongly coupled to pn by p exchange (see Figure [[]). So 
we consider the process p + —> b\n —> pn + . 


For on-shell states the Lorentz invariant amplitudes can be parameterized as 



f pM* 

C /i C ^ 


(9p bl n9^ + g^Pb^p) 


( 1 ) 


M 


ftl—>P' 


_ M p* 
'V t p 


(9b 1P y v + g^pvPpPbJ 


( 2 ) 


Mp-n+„o = g(m^{jpP + - P%) (3) 

where Px and e * refer to the momentum and polarization 4-vectors of X respectively, 
is the flat space metric tensor, and g% bl7V , gf bin , g bl ™, g bipr] and g pn7r are decay constants to 
be determined. These are discussed in Appendix B. 

The doorway amplitude for the process p — » b^n 0 gn + is then 
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M 


doorway 


7^ / ^Pbl e i^ S pbl^ U + g^nPh^p)— - 9UtT 2 I ■ 

(27r) 4 J ^ H Pbi~ m b\ + ze 

1 


X istm^ + 9b lPV PpPb ! 


“2 ^2 i - “ P^o) - j- ^ . (4) 

Pp - + «e p^ 0 - m^o + 76 


where e is a small real number. Here we have contracted the Lorentz indices on the internal 
vector particles in the usual way m , effectively working in the “Feynman gauge”. Inte¬ 
gration is performed over the loop momentum. We evaluate the doorway amplitude in the 
rest frame of p using conservation of momentum at the vertices and the Cutkovsky rules to 
obtain [|H| 


Mr 


^gp-n 


{(si + S 2 £ + S 3 z 2 ) In 


1 + z, 


2(s 2 + ZS 3 )} 


(5) 


doorway 32n(E bl + E^)q 

if the component of the polarization of p is in the direction of the outgoing particles, i.e. p 
or 7r + . For other polarizations, Mdoorway = 0. 

Here 


Si — —gpbmQ (Pbipp + gblppiEbiiEnO + E n +) + p 2 )} 


( 6 ) 


s 2 = p {-g^g^pv + g^g^m^E^ + K+)+ 

(g^gSpr, + g? bl „g°p V mp{E bl - E v ))(E bl (E n o + E n+ ) + p 2 ) - g% in g bipri q 2 } (7) 


s 3 = gb 1 PV p 2 q {g S p bl n + gS> lv m p ( E bl - e v )} 

p 2 + q 2 - (E b 1 - E v ) 2 + m 2 

z =-- 

2 pq 


( 8 ) 

( 9 ) 


and p is the magnitude of the b\ or 7T° momentum, q the magnitude of the p or 7 r + momen¬ 
tum, and Ex the energy of X, all in the rest frame of p. 

The doorway decay width for p + —> bin —> pn + is then 












r 


P — >b\TT —>7]7T+ 


87T?7l? 3 


12 M 


doorway \ 


(10) 


where we have taken into account that there are two possible intermediate processes con¬ 
tributing to the total amplitude, i.e. p + —> bfir 0 — > pir + and p + —> —> ?/7r + . 

We calculate that the doorway width is 57 ± 14 MeV for all the particles on-shell, taking 
into account uncertainties in the couplings bi — > p + p and p + — * 7T + 7r° (see Appendix |B]). 
As remarked in Appendix |B], there are uncertainties in the coupling p + —*■ bf 7T° which can 
make this doorway width up to ~ 40% smaller. Thus we conclude that the doorway width 
is less than 57 ± 14 MeV. It should also be remembered that the doorway calculation as it 
stands provides an upper limit, since we would get a smaller answer if we were to take one 
of the internal legs off-shell and introduce form factors |^0|. However these are unknown, 
and as 57 ± 14 MeV is well above the pn width required this does not create a problem. 


3.2 Non—resonant ijtt Deck background 

The 1.4 GeV peak in the pit channel occurs in the vicinity of the pn and bin thresholds, 
and it is therefore natural to consider these as being responsible in some way for the pn 
peak. The Deck mechanism |19| is known to produce broad low-mass enhancements for a 
particle pair in three-particle final states, for example in np —> (pn)p. In this latter case, 


the incident pion dissociates into pi r, either of which can then scatter off the proton 


At sufficiently high energy and presumed dominance of the exchange of vacuum quantum 
numbers (pomeron exchange) for this scattering one obtains the “natural parity change” 
sequence n —> 0 - ,l + ,2 - .... (the Gribov-Morrison rule |22|]). However if the scattering 
involves the exchange of other quantum numbers then additional spin-parity combinations 
can be obtained, including J p = 1 _ . This can be seen explicitly in ref. [H| for the 
reaction np —> (pn)p in which the full up scattering amplitude was used, so that the effect 
of exchanges other than the pomeron are automatically included. The J p sequence from 
the “natural parity change” dominates due to the dominant contribution from pomeron 
exchange, but other spin-parity states are present at a non-negligible level. The Reggeised 
Deck effect can simulate resonances, both in terms of the mass distribution and the phase 


23|j. It can produce circles in the Argand plot, the origin of which is the Regge phase 


factor exp[— i^na(tn)]. 
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Figure 2: Deck background production in pn. 
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It is also important to note that rescattering of the lighter particle from the dissociation of 
the incident beam particle is not a prerequisite, and indeed both can contribute |2T|. We 
suggest that in our particular case the relevant processes are (from left to right in Figure 

® 


1.7r —> biu, up —> up giving a bin final state. 

2. 7T —> 7rp, pp — »■ rjp giving a pn final state. 

3. 7T —> p7r, 7rp —> 7rp and pp —> pp giving a Tip final state. 

For each of these processes the rescattering will be predominantly via p (natural parity) 
exchange to give the required parity in the final state. Obviously process (ii) produces 
a final pn state directly, but for (i) and (iii) the &i7r and p7r final states are required to 
rescatter into pn (for which the doorway calculation provides an explicit mechanism). 

Unfortunately only the up cross section can be obtained with any reliablity. The others 
can be estimated with varying degrees of uncertainty from: 

1. up —> b\ tt: data on 7r°p —> un, which can be inverted to give up —> 7r°p. 

2. Tip — > pn: data on pp — > pp , which by assuming vector meson dominance can give 

pp —> pp. 

3. Tip — > pn: data on Ti~p — > Ti°n and on 7r =b p — >• 7r ± p; data on pp —> pp and vector 
meson dominance. 

In view of the uncertainties in the underlying reactions, the lack of an explicit value for 
the J PC = 1“ ' + pn cross section in the E852 experiment, and the impossibility of a precise 
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evaluation of the rescattering into the pn channel from pn and bin, we have not attempted 
a complete Deck-type calculation. We concentrate rather on the mass-dependence which 
it generates. The characteristic mass-dependence is a peak just above the threshold. Thus 
there are three peaks from our proposed mechanism: a sharp peak just above the pn thresh¬ 
old; a broader one at about 1.2 GeV from the pn channel; and a very broad one at about 
1.4 GeV from the b\n channel. The first of these is effectively removed by experimental 
cuts, but the net effect of the two latter is to produce a broad peak in the pn channel. 
Thus invoking this mechanism does provide an explanation of the larger width of the pn 
peak at 1.4 GeV in the E852 data compared to that of the pn peak at 1.6 GeV. Because 
of the resonance-like nature of Deck amplitudes it is also possible in principle to simulate 
the phase variation observed. However as there are Deck amplitudes and the 1.6 GeV 
resonance, presumably produced directly, it is necessary to allow for interference between 
them. We use the K-matrix formalism to calculate this, and also to demonstrate that the 
Deck mechanism is essential to produce the 1.4 GeV peak. 


3.3 K—matrix with P—vector formalism 


It is straightforward to demonstrate that within the K-matrix formalism it is impossible 
to understand the pn peak at 1.4 GeV as due to a 1.6 GeV state if only resonant decays to 
pn, pn and bin are allowed despite the strong threshold effects in the two latter channels 
[]. We find that for a bin width of ~ 200 MeV and pn and pn widths in the region 1 — 200 
MeV there is no shift of the peak. However, when a non-resonant pn P-wave is introduced, 
the interference between this and the 1.6 GeV state can appear as a 1.4 GeV peak in pn. 


We have seen that the non-resonant pn wave can have significant presence at the bin or fin 
threshold (called the “P+S” threshold), e.g. 1.368 GeV for bin, because of the substantial 
“width 1 ' generated by the Deck mechanism. Since the hybrid is believed to couple strongly 


to “P+S” states due to selection rules [JT3|, [15[], the interference effectively shifts the peak 
in pn down from 1.6 GeV to 1.4 GeV. It is not necessary for the 1.6 GeV resonance to 
have a strong pn decay. It is significant that the E852 experiment finds p at 1370 ± 16i|o 
MeV, near the bin threshold, but not at 1.6 GeV. It is possible for a state to peak near 
the threshold of the channel to which it has a strong coupling, assuming that the (weak) 
3 The use of bin is not critical here: any channel with a threshold near 1.4 GeV will suffice. 
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channel in which it is observed has a significant non-resonant origin. 


We follow the K-matrix formalism in the P-vector approach as outlined in | [24] . |25| . We 
assume there to be a p with rrip = 1.6 GeV as motivated by the structure observed in pi r 
[3j. The problem is simplified to the case where there is decay to two observed channels 
i.e pn and /nr, and one unobserved P + S channel. These channels are denoted 1, 2 and 
3 respectively. The production amplitudes and the amplitude after final-state interactions 
are grouped together in the 3-dimensional P- and F-vectors respectively. In order to 
preserve unitarity [Q we assume a real and symmetric 3x3 K-matrix. The amplitudes 
after final-state interactions and production are related by |24 


F = (/ — iK) P 


( 11 ) 


We define the widths as 


r* 



B\gi) 


p(Qi) 


i = 1,2 


( 12 ) 


r 3 = 73 r p p(Qa) (I 3 ) 

where g* is the breakup momentum in channel i from a state of effective mass w, and gf is 
the breakup momentum in channel i from a state of effective mass mp. The kinematics is 
taken care of by use of the phase space factor 


2 q 

w 

and the P-wave angular momentum barrier factor 



(14) 


B\q) 


0 q/q R ) 2 

i +( q/qR ) 2 


where the range of the interaction is qn = 1 fm = 0.1973 GeV. 


(15) 


We assume the experimental width in pi: of Yp = 168 MeV |3[| to be the total width of 
the stateQ. We adopt the flux-tube model of Isgur and Paton and use the pi r and biir 
4 It is found that our results in Fig. 3 are very similar even for a width of 250 MeV. 
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widths which it predicts for a hybrid of mass 1.6 GeV. Since the model predicts that the 
branching ratio of a hybrid to biir is 59 — 74 % and to /i7r is 12 — 16 % [17], we obtain the 


P + S'-wave width to be 120 — 150 MeV. Analysis of the data shows that the /nr branching 
ratio of p(1600) is 20 ± 2 % ||, corresponding to a pi r width of 30 — 37 MeV. This is 
consistent with flux-tube model predictions of 9 — 22 % |L7| . For the simulation we use a 


&i7r width of 120 MeV, a pir width of 34 MeV, and an pir width of 14 MeV, well within the 
limits set by the doorway calculation. We neglect other predicted modes of decay since we 
restrict our analysis to three channels. 


The K-matrix elements are 


Kij = 


m 2 ~ — w 2 


+ c. 


V 


(16) 


where c tJ includes the possibility of an unknown background. 

In the simulation we assume that the Deck terms can be treated as conventional resonances. 
This is not necessary, but is done to reduce the number of free parameters. We assume 
that the r/7r Deck amplitude is produced predominantly via the bpn and pir channels, and 
so is modelled as a resonance at a mass mb i = 1-32 GeV and a width T&i = 300 MeV. This 
width fits the E852 data at low pn invariant masses (see Figure [|a). The pir background 
is assumed to peak at a mass nib 2 = 1-23 GeV with a width r& 2 = 400 MeV, which when 
plotted as an invariant mass distribution effectively peaks at ~ 1.15 GeV, in agreement 
with detailed Deck calculations in the 1 ++ wave 0- 

We incorporate the pn and pir Deck background by putting c t] = 0 except for 


_ ITlbl^bl _ m b2^b2 

Cn — —5 - 5 - C 22 — —5 -; 

— w z 2 — w 


The widths are defined analogously to Eq. [L2] as 


(17) 


I’m = Tm r p > = 1.2 (18) 

where q\ is the breakup momentum from a state of effective mass (for i — 1, 2). 

The production amplitudes are given by 
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where the (dimensionless) complex number Vp measures the strength of the production of 
p. We take C 3 = 0 and 


Cl = 


m bl v bly Jr bl ri 

m'l x — w 2 


c 2 = 


rrib2Vb2\J^b2F f, 

m b2 ~ W2 


( 20 ) 


where the complex numbers V&j gives the production strengths of the Deck background in 
channel 1 . 


The results of this fit are shown in Fig. |3] and clearly provide a good description of the pn 
data |L], |25 . 


We briefly discuss the results. Fig. |a indicates a steep rise for low invariant pir masses, 
and a slow fall for large pn masses. This naturally occurs because of the presence of the 
resonance at 1.6 GeV in the high mass region, which shows as a shoulder in our fit. Figure 
|3]b reproduces the experimental slope and phase change in pn |25|]. One might find this 
unsurprising, since the background changes phase like a resonance. However, we have 
confirmed, by assuming a background that has constant phase as a function of pn invariant 
mass, that the experimental phase shift is still reproduced. The experimental phase shift 
is hence induced by the resonance at 1.6 GeV. 


Our fit to E852 ip r and pi r data (with a prediction for the bin data) requires 12 independent 
parameters (see the caption of Figure 3). 

Without the inclusion of a dominant P + S'-wave channel [] the pn event shape clearly shows 
two peaks, one at 1.3 GeV and one at 1.6 GeV, which is not consistent with the data [ 1 ]. 
The phase motion is also more pronounced in the region between the two peaks than that 
suggested by the data ^5[. The role of the dominant P + A-channel is thus that at invariant 
masses between the two peaks, the formalism allows coupling of the strong P + S channel 
to pn, so that the pn appears stronger than it would otherwise, interpolating between the 
peaks at 1.3 and 1.6 GeV, consistent with the data HJ. A dominant P + S decay of the p 
is hence suggested by the data. 


5 The &i7r coupling of the resonance is set to zero, with the rjir and pn couplings the same as before. 

















4 Discussion 


We have argued that on the basis of our current understanding of meson masses it is 
implausible to interpret the 1.4 GeV peak seen in the J PC = 1 _+ pir channel by the BNL 
E852 experiment as evidence for an exotic resonance at that mass. We acknowledge that 
this is not a proof of non-existence and note the Crystal Barrel claim for the presence of 
a similar state at 1400 ± 20 ± 20 MeV in the reaction pp —> pir + n~ . However this is not 
seen as a peak and is inferred from the interference pattern on the Dalitz plot. It has 
not been observed in other channels in pp annihilation at this mass, which is required for 
confirmation. So at present we believe that the balance of probability is that the structure 
does not reflect a real resonance. 

Given this view, it is then necessary to explain the data and in particular the clear peak and 
phase variation seen by the E852 experiment. Additionally the observation of the peak only 
in the pn channel, which is severely suppressed by symmetrization selection rules, requires 
justification. We have dealt with these two questions in reverse order. We first demonstrate 
that final-state interactions can generate a sizable pi r decay. We believe that this result 
by itself is of considerable significance and is of wider relevance. We then suggest that the 
E852 pn peak is due to the interference of a Deck-type background with a hybrid resonance 
of higher mass, for which the p at 1.6 GeV is an obvious candidate. This mechanism also 
provides the natural parity exchange for the former which is observed experimentally. The 
parametrization of the Deck background is found not to be critical. 

A key feature in our scenario is the presence of the large “P + S v amplitude which drives 
the mechanism. This should be observable both as a decay of the 1.6 GeV state and as a 
lower-mass enhancement due to the Deck mechanism. Depending on the relative strength 
of these two terms the resulting mass distribution could be considerably distorted from a 
conventional Breit-Wigner shape as the Deck peak is broad and the interference could be 
appreciably greater than in the pn channel. 
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A Appendix A: Decay of orbitally excited hybrid 


We detail here the flux-tube model calculation for the decay of an orbitally excited hybrid 
to P + S -wave mesons. The normalized wave functions of the P-wave and S-wave mesons 
are just S.H.O. wave functions with the same inverse radius f3 |U|, and are respectively 


2 /f ^174 rY iM L d) e xp ~/3 2 r 2 


^ 1,2 2 

“V 4 exp--/3 r 


The normalized wave function of the orbitally excited hybrid is taken to be 


( 21 ) 


Mp r&vL M* a( 0> “0) ex P ~2^ r2 




^ = 


(2 L + l)/3| +2<5 
\ 27rr(| + h) 


( 22 ) 


where the Wigner D-function guarantees that the state has total orbital angular 

momentum L — 2, the first orbital excitation above the ground state with total orbital 
angular momentum 1, is the total orbital angular momentum projection, and A = ±1 
is the angular momentum of the flux-tube around the QO axis 0. The inverse radius 
f3p characterizes the size of the wave function and T is the Gamma-function. The radial 
dependence is chosen to be proportional to r s , where 5 is chosen such radial Schrodinger 
equation [|TG|] in the limit r —> 0, which leads to the condition 6(6 + 1) = L(L + 1) — A 2 [jl3| , 
implying that 6 = 1.79 for L — 2. The lowest orbitally excited hybrid has the QQ in spin 
1, just like the ground state hybrid. 

The relevant overlap can be obtained by inserting the spacial wave functions into the de¬ 
cay matrix element and performing the integration over the quark-antiquark pair creation 


position |L3 


0 62 yn \/2 /*2-n‘ /*7r roo 

M m£m l = flavour ( 1+ alp i d<i> J 0 s[ll9d9 r2 dr r5 + P) ' 

x v Mi i(^> -0) v1 m l i(^> -0) ex P ^p • r ex P -y r2 ( 2 ^I + P 2 ) ( 23 ) 

where p is the momentum of the outgoing mesons in the rest frame of the hybrid and e 
a spherical basis vector. Notice that the pair creation constant 70 of the 3 Pq model enters 


16 






















explicitly in Eq. |23|. This is because the flux-tube model, within the assumptions made for 
the wave functions, gives a prediction for the couplings of a hybrid in terms of couplings 
for mesons in the 3 P 0 model |13| (the constants 0.62 and 0.1 in Eq. 6 are derived from 


flux-tube dynamics). We take y 0 = 0.39 Jl3[. ^6|, [HJ]. The integral in Eq. 
numerically. 


is performed 


We can write the decay amplitudes in terms of the amplitudes in Eq. ^3| as follows, following 
ref. 0. For b x i r (“flavour” = 2): 


S-wave amplitude =- 7 =^ (76A1 2 i + 73Adn + -Aden + 73A1io + Aloo) 


D-wave amplitude = 


VlE 

1 

730 


A (y/6M 2 i + y/3M n + M 01 - 2V3M 10 - 2 M 


00 


where A selects the imaginary part of the amplitude. For fill and airj (“flavour” = 
1 respectively): 


) (24) 

72 and 


S-wave amplitude 
D-wave amplitude 


—^=9 (76Ad 21 + 73Ad 11 + Ad 0 i + 73A1 io + .Ad 00 ) 

^==9 (76Af 2 i — 273 Ad n — 5Af 0 i T 73Afio + Afoo) (25) 


The A'i(1270) is regarded as cos 6 | J Pi) + sin 9 | 3 Pi) and Ad(1400) the orthogonal partner, 
with 9 = —34° |T3]. |T(J. 1 Pi and 3 Pi are the P-wave mesons with QQ combinations of the 
decay amplitudes to the 1 P i meson (Eq. ^4 ]) and 3 Pi meson (Eq. ^5] ). For Ad (1270)A' and 
Ad (1400)A" “flavour” = 72. 

For / 2 7 t, a 2 7? and ATI(1430 )K (“flavour” = 72, 1 and 72 respectively): 


D-wave amplitude = —(76Af 2 i — Aloi — VSM .10 — Afoo) (26) 

275 


The decay amplitudes in Eqs. 


are then used to calculate widths according to the 


phase space conventions of Eq. 6 of ref. [13 


B Appendix B: Doorway calculation constants 
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B.l 


9 S phn and 9 %r 


These couplings cannot be obtained from experiment, as there is currently no published 
data on the p + coupling to bf n°. We have calculated g~ bin and gf bin in the non-relativistic 


flux-tube model of Isgur and Paton, following the conventions and methods of ref. |13 


(except that we assume the relativistic phase space convention |26|]). The wave functions 
of the mesons and the hybrid are given respectively by Eqns. ^T] and [22] (with L — 1). We 
find that 


( s 

ypb\7 
\ 9pb\7 


where 


\® 3 / 




= -J im pE bl E n o 


m b 
\ ™.pP 


— 2zui + 3vj2 + ^3 

^(3 {vj -2 + zu 3 ) + (^— 1)(— 2m 1 + 3zu 2 + G 7 3 )) 


/ 


O- 62 70 2 n Pp + f°° dr r 2+5 cxpf r ‘ (28 2 I 8 2 )) 

(l + ^) 2 /5^3r(| + h)io drr 6XP( 4 ( ^ +P)) 


P 2 rj 0 ( p i) 


\ 


P3 i(f) 

V /? 2 j 2 (f) j 


(27) 


(28) 


where j t and T refers to the spherical Bessel and Gamma functions respectively. We use 
70 = 0.53 which reproduces conventional meson decay phenomenology for relativistic phase 
In Eq. |27|, (3 refers to the inverse radius of the bf or 7T°, the parameter that 


space 


26 


enters in the S.H.O. wave function. Similarly, (3p is the inverse radius of p + . 
Setting S — 1, (3p — 0.27 GeV and [3 = 0.36 GeV, yields 


= 3.0 GeV = -8,2 GeV 


7—1 


(29) 


We chose a value of y 0 towards the upper end of the range in the literature. In calculations 
of excited mesons, values of y 0 as low as .4 have been used P7| . The values for g^ bin and 
gj? bin can hence be only .4/.53 of the values quoted. 


B.2 




9b\pp aad 9h\ prj 


Although the iq coupling to p + g is not known experimentally [g 8 J, its coupling to uir + 
is well known [^jj, and can be used to obtain the p + r/ coupling. We first derive the urn 
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coupling by assuming that 100% of the decays of b± is to u>tt + and using the experimentally 
measured D-wave to S-wave amplitude ratio. 

The amplitude for —> u;7r + can be written as 

M bl ^ w = + g^PuPbJ ( 30 ) 

Using the Jacob-Wick formulae we write the S-wave and D-wave decay amplitudes as 


M s bl _ 



Vz { 


-zgL 


E u _ 


m u 


dbiujiz 


+ Pl—9b-_ U) ( 31 ) 


m„ 


= yf ((E - l)gU - (32) 

where all energies and momenta refer to the bf rest frame. Relating the S-wave and D- 
wave amplitudes to the corresponding widths in the usual way (analogous to Eq. 0 , we 
finally obtain 


m bl (V 2 + 





(33) 


Using the experimental data T tota i — 142 ±8 GeV and y y— — +0.29 ±0.04 (where the sign 
is taken to mean that the D-wave and S-wave amplitudes have the same sign), we obtain 
can then deduce the coupling constants 


9b,um — 4.6 ±0.2 GeV 9 ,% = 14.4 ±2.2 GeV- 1 (34) 

To obtain the b\ coupling to p + rj, we note that (neglecting effects due to phase space), 
it should be related to the u;7r + coupling by a simple flavour factor. This is because the 
Lorentz structure of the two decays are identical. Assuming that the decomposition of the 
g which is motivated by experiment, i.e. g = ^(^(wh + dd) + ss ), we have 

s;„ = 44-. = 3 - 2 ±0.1GeV 9b,„=^=9t,= W.2 ±1.5 GeV-' (35) 
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We have also performed a flux-tube ( 3 P 0 ) model calculation to independently derive the 
coupling constants. For y 0 = 0.53, /3 = 0.4 GeV we obtain 


= 3.4 GeV 


< m = 10.6 GeV- 1 


(36) 


The agreement (both in sign and magnitude) is clearly impressive, underlining the signifi¬ 
cant agreement of the 3 P 0 model with experiment [^6] . 


B-3 (jp T 


Here we assume that 100% of the decays of p + are to 7r + 7r° |28|. We evaluate the amplitude 
Eq. in the rest frame of p + and connect the amplitude to the width of 150.7 ± 0.6 MeV 
HI , according to the usual relation (analogous to Eq. m to obtain 


<W = 6.02 ± 0.02 


(37) 


References 

[1] D.R. Thompson et al. (E852 Collab.), Phys. Rev. Lett. 79 (1997) 1630. 

[2] A. Abele et al. (Crystal Barrel Collab.), Phys. Lett. B423 (1998) 175; W. Diinnweber 
(Crystal Barrel Collab.), Proc. of HADRON’97 (Upton, N.Y., August 1997), p. 309, 
eds. S.-U. Chung, H.J. Willutzki. 

[3] N.M. Cason (E852 Collab.), Proc. of Intersections between Nuclear and Particle Physics, 
p. 471 (Big Sky, Montana, May 1997), ed. T.W. Donnelly, American Institute of Physics; 
D.P. Weygand and A.I. Ostrovidov (E852 Collab.), Proc. of HADRON ’97 (Upton, N.Y., 
August 1997), p. 313, 263, eds. S.-U. Chung, H.J. Willutzki; A.I. Ostrovidov et al. (E852 
Collab.), in preparation. 

[4] Yu.P. Gouz. et al. (VES Collab.), Proc. of the 26 th ICHEP (Dallas, 1992), p. 572, ed. 
J.R. Sanford. 

[5] D. Ryabchikov (E852 Collab.), “Study of the reaction ir~p —> p in E852 exper¬ 

iment”, Proc. of HADRON'97 (Upton, N.Y., August 1997), p. 527, eds. S.-U. Chung, 


20 





H.J. Willutzki; Yu.P. Gouz. et al. (VES Collab.), Proc. of the 2Q th ICHEP (Dallas, 
1992), p. 572, ed. J.R. Sanford. 

[6] P.R. Page, Phys. Lett. B415 (1997) 205. 

[7] T. Barnes, F.E. Close and E.S. Swanson, Phys. Rev. D52 (1995) 5242. 

[8] C. Michael et al ., Nucl. Phys. B347 (1990) 854; Phys. Lett. B142 (1984) 291; Phys. 
Lett. B129 (1983) 351; Liverpool Univ. report LTH-286 (1992), Proc. of the Workshop 
on QCD : 20 Years Later (Aachen, 9-13 June 1992). 

[9] P.R. Page, Phys. Lett. B401 (1997) 313; F. Iddir et al., Phys. Lett. B207 (1988) 325; 
H.J. Lipkin, Phys. Lett. 219 (1989) 99; S.F. Tuan, T. Ferbel, R.H. Dalitz, Phys. Lett. 
B213 (1988) 537; S.F. Tuan, Mod. Phys. Lett. A, Vol.6, No.9 (1991) 789. 

[10] N. Isgur and J.E. Paton, Phys. Rev. D31 (1985) 2910; R. Kokoski, N. Isgur, Phys. 
Rev. D35 (1987) 907. 

[11] S. Collins (UKQCD Collab.), Nucl. Phys. Proc. Suppl. 63 (1998) 335. 

[12] 1.1. Balitsky, D.I. Dyakonov, A.V. Yung, Z. Phys. C33 (1986) 265. 

[13] F.E. Close, P.R. Page, Nucl. Phys. B443 (1995) 233; Phys. Rev. D52 (1995) 1706. 

[14] F. Iddir, S. Sahr, O. Pene, |icp-ph/9803470 . 

[15] P.R. Page, Phys. Lett. B402 (1997) 183. 

[16] J.I. Latorre et al., Z. Phys. C34 (1987) 347; S. Narison, “QCD spectral sum rules”, 
Lecture Notes in Physics, Vol. 26 (1989), p. 369; J. Govaerts, F. de Viron, Phys. Rev. 
Lett. 53 (1984) 2207. 

[17] P.R. Page, Proc. of HADRON’97 (Upton, N.Y., August 1997), p. 121, eds. S.-U. Chung, 
H.J. Willutzki, |hcp-ph/9711241| ; P.R. Page, E.S. Swanson, A.P. Szczepaniak, in prepa¬ 
ration. 

[18] M.P. Locher, Y. Lu and B.S. Zou, Z.Phys. A347 (1994) 281. 

[19] G. Ascoli et al. Phys. Rev. D8 (1973) 3894; ibid. D9 (1974) 1963. 

[20] O. Gortchakov et al., Z.Phys. A353 (1996) 447. 


21 







[21] L. Stodolsky, Phys. Rev. Lett. 18 (1967) 973. 

[22] V. Gribov, Sov. J. Nucl. Phys. 5 (1967) 138; M.R.O. Morrison, Phys. Rev. 165 (1968) 
1699. 

[23] M.G. Bowler et al ., Nucl. Phys. B97 (1975) 227; M.G. Bowler, J. Phys. G5 (1979) 
203. 

[24] I.J.R. Aitchison, Nucl. Phys. A189 (1972) 417; S.-U. Chung et al., Ann. der Phys. 4 
(1995) 404. 

[25] S.-U. Chung, Private communication; Proc. of Jefferson Lab / NCSU Workshop on 
Hybrids and Photoproduction Physics (1997). 

[26] P. Geiger, E.S. Swanson, Phys. Rev. D50 (1994) 6855. 

[27] T. Barnes, F.E. Close, P.R. Page, E.S. Swanson, Phys. Rev. D55 (1997) 4157. 

[28] Particle Data Group, Phys. Rev. D54 (1996) 1. 


22 



Figure 3: Results of the K-matrix analysis, (a) The events (|F\| 2 ) in r/n as compared to 
experiment [|l|; (b) The phase (of Fi) in r/n compared to experiment |25|]. The invariant 


mass w is plotted on the horisontal axis in GeV. When the phase is plotted it is in radians, 
with the overall phase ad hoc. The parameters of the simulation are trip = 1.6 GeV, 
Tp = 168 MeV [|], 71 = 0.31, 72 = 0.52, 73 = 1.49, m bl = 1.32 GeV, m b2 = 1.23 
GeV, 7 bi = 1.53, 7 ^ = 2.02, V b \/Vp = 2.05e 2 ' 77 *, V^/Hi = 0.35e 16 *. Vp sets the overall 
magnitude and phase, which is not shown. None of the ratios of production strengths 
should be regarded as physically significant, since the K-matrix formalism allows for the 
introduction of additional parameters in the modelling of the backgrounds, which would 
change the values of these ratios. The plots shown here are only weakly dependent on the 
pir parameters 7^2 and V b 2 - The parameters have been chosen to fit both the pn data |]| 
and the preliminary /wr data [31]. Experiment has not been able to eliminate the possibility 
that the low mass peak in pi r is due to leakage from the a\. The background amplitude in 
pir is being used as a means of parametrising all forms of background into the pir channel, 
including leakage or Deck. 
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